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Initial results on particle yields obtained by the HARP experiment are presented. The measurements correspond 
to proton-nucleus collisions at beam energies of 12.9 GeV/c and for a thin Al target of 5% interacion legth. The 
angular range considered is between 10 and 250 mrad. This results are the first step in the upcoming measurement 
of the forward production cross-section for the same target and beam energy, relevant for the calculation of the 
far-to-near ratio of the K2K experiment. 



1. Introduction 

The HARP experiment \, was designed to per- 
form a systematic and precise study of hadron 
production for beam momenta between 1.5 and 
15 GeV/c and target nuclei ranging from hydro- 
gen to lead. The detector was located at CERN, 
in the PS beam.The DAQ recorded 420 million 
events during the years 2001 and 2002. 

The physics program of HARP includes: a) the 
measurement of pion yields for a variety of ener- 
gies and targets relevant for he design of the pro- 
ton driver of a future neutrino factory b) the 
measurement of pion yields on low Z targets as 
well as on cryogenic oxygen and helium targets, 
useful to improve the precision of atmospheric 
neutrino flux calculations 0; and c) the measure- 
ment of pion and kaon yields, relevant for the 
calculation of the neutrino fluxes of experiments 
such as MiniBooNEg] and K2K0. 

HARP (Fig. Pi is a large acceptance spectrom- 
eter, with two distinct regions. In the forward 
part of the apparatus (up to polar angles of about 
250 mrad), the main tracking devices are a set 
of large drift chambers. Magnetic analysis is pro- 
vided by a 0.4 T dipole magnet and particle iden- 
tification relies in the combination of a threshold 
Cherenkov detector, a time-of- flight wall and an 
electromagnetic calorimeter. In the rest of the 



solid angle the main tracking device is a TPC, 
which is complemented by a set of RPC detectors 
for time-of-flight measurements. The target is lo- 
cated inside the TPC. In addition, sophisticated 
beam instrumentation (including three timing de- 
tectors and threshold Cherenkov detectors) pro- 
vides identification of the incoming particle and 
allows the interaction time at the target to be 
measured. 




Figure 1. Schematic layout of the HARP spec- 
trometer. 
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Given the immediate interest of the Mini- 
BooNE and K2K experiments in a measurement 
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of the production cross sections for pions and 
kaons at the energies and targets relevant for their 
beam setups, the HARP collaboration has given 
priority to the analysis of those particular data 
sets. In particular, we present in this article an 
initial result relevant for the K2K experiment. 
Specifically, we measure pion yields from proton- 
nucleus collisions, for a thin Al target ( 5%A ) and 
a beam energy of 12.9 GeV. The K2K neutrino 
beam is produced from the decay of hadrons (pri- 
marily pions) emanating from a 2A Al target, at 
the KEK PS proton driver energies of 12.9 GeV. 
Further details about the motivation of this anal- 
ysis can be found in Ref. [§J. 

2. Forward Tracking and Particle Identifi- 
cation 

2.1. Tracking 

Tracking of forward-going particles is done by 
a set of large drift chambers(NDC) placed up- 
stream and downstream of the dipole magnet. 
The chambers were recuperated from the NO- 
MAD experiment and their properties have been 
described elsewhere[7]. Each NDC module is 
made of four chambers, and each chamber of three 
planes of wires with tilted angles —5°, 0° and 5°. 
The single-wire efficiency is of the order of 80%, 
and the spatial resolution approximately 340 \ira. 
The performance of the NDC is inferior to the one 
measured in NOMAD, where the hit efficiency 
was close to 95% and the resolution approached 
150 fim. This is mainly due to the fact that, for 
HARP, the NDC used a gas mixture and voltage 
settings different from the ones in NOMAD 2 . 

Due to the low hit efficiency, the algorithm used 
in NOMAD, based in building up space points 
(with an efficiency of about .95 3 = 85%) is not 
viable (for Harp we would obtain .80 3 = 51%). 
Instead, the reconstruction builds 2D and 3D 
track segments in each NDC module (12 hits max- 
imum), which are fitted to a straight line model 
via a Kalman Filter fit|U]. 

Next, the algorithm attempts all possible com- 
binations (which include at least a 3D segment) 
to connect tracking objects in the modules down- 

2 Ar(90%)-CO 2 (9%)-CH 4 (l%), a non-flamable gas mix- 
ture. 



stream of the dipole magnet. Combinations such 
as 3D + 2D or even 3D + hits-not-associated 
are valid ways to build longer, 3D segments. 

To measure the momentum it is necessary to 
connect a (3D) segment downstream of the dipole 
with at least one space point upstream of the 
dipole. Since one can impose the constraint that 
all tracks emanate from the event vertex 3 this 
point is always known and therefore the necessary 
and sufficient condition for a particle emanating 
from the target to have its momentum measured 
is that a 3D segment can be measured by the 
combination of downstream NDC chambers. 

A measurement in the first NDC module is not 
strictly necessary to analyse the track momen- 
tum. Ideally, one would like, of course, to connect 
a 3D segment downstream with a 3D segment in 
the NDC1 and the vertex point. However, the 
tracking efficiency of the upstream module is size- 
able lower than the efficiency of the downstream 
modules. This is due to, a) a higher hit density in 
NDC1 4 and b) the lack of redundance due to the 
fact that there is only a single module upstream of 
the dipole. A too restrictive condition to accept 
a track such as matching a 3D segment upstream 
of the dipole with a 3D segment downstream re- 
sults not only in lower efficiency, but also in a 
systematic error, arising from the uncertainties 
in the hit density expected in NDC1 which would 
be hard to estimate (hit density depends on fac- 
tors such as track multiplicity and opening angle, 
which are model-dependent). Instead, by build- 
ing a track always when a 3D segment is present 
downstream of the dipole, one is practically inde- 
pendent of the tracking efficiency in NDC1. 

In order to quantify these considerations, we 
have considered three different types of tracks, 
depending on the matching between the down- 
stream and upstream modules. Type I is a 3D-3D 
matching. Type II matches a 3D (downstream) 
with a 2D (upstream). Type III are those where 
one has been unable to find either a 3D or a 2D 



3 for a thin target the vertex is known with a precision of 
about 1 cm In the case of thick target, a vertex algorithm 
will provide the event vertex to a comparable precision. 
4 the hit density decreases quadratically with the distance 
to the vertex, and therefore the effect is only significant 
for NDC1. 
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segment in NDC1 and thus the track is built us- 
ing the downstream modules and the vertex con- 
straint. 

2.2. Particle Identification 

Particle identification (PID) in the forward 
region of the spectrometer combines the infor- 
mation provided by beam detectors and three 
systems located downstream the dipole mag- 
net. Namely, a threshold Cherenkov detector 
(CKOV), a time-of-fiight wall (TOF), and an 
electromagnetic calorimeter (electron identifier, 
EID). These subsystems have been described in 
Ref. [H]. Its combined information results in 
good PID over the whole range of relevant mo- 
menta, as well as redundancy due to overlaps. 
Pion/proton separation is provided by TOF up 
to 4.5 GeV/c, and by the CKOV above 3 GeV/c. 
Electron/pion separation is covered by the CKOV 
below 3 GeV/c and by the EID above 2 GeV/c. 
Finally the kaon contamination can be estimated 
with the CKOV above 3 GeV/c and with the 
(TOF) below this energy. 

3. The analysis 

About one sixth (1 million events) of the "K2K 
thin target" data has been analysed. The unnor- 
malised pion production differential cross section 
can be computed as follows: 
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where the sum runs over track types (I = 3D- 
3D, 11= 3D-2D, III=3Z?-vertex) and the indices 
i and j correspond to true and reconstructed 
(p, 8) bins respectively. ef cc is the geometrical 
acceptance, e t j rack is the tracking efficiency, Mij 
is the migration matrix from reconstructed bin 
j to true bin i, ej is the pion identification effi- 
ciency, t]J is the pion purity and N? is the ob- 
served pion yield. The pion purity is defined as 
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V /NJ = (N7 - N° K9 )/NJ, where 
j^true-n j g number of true observed pions in 
the bin j, and N^ kg is the number of particles 
misidentified as pions in the same bin j. That 
means that 77J is just the probability of correct 



identification of a pion. The pion identification 
efficiency, ej is the fraction of times that a tracked 
pion is identified as such. 

Fig. |21 shows the spectrometer acceptance. For 
clarity of ilustration, the momentum acceptance 
has been computed in the region of good angular 
acceptance, {\9 y \ < 50 mrad, \9 X \ < 200 mrad) 
and conversely, the geometrical acceptance has 
been computed in the region of good momen- 
tum acceptance (P > lGeV). As it can be seen 
the acceptance is good for pion momenta above 1 
GeV/c, and still acceptable above 0.5 GeV. As 
illustrated in Ref. [H] this permits the measure- 
ment of K2K neutrino fluxes up to 0.25 GeV, 
which amply covers the energy region where the 
atmospheric oscillation affects the neutrino spec- 
trum maximally. The angular acceptance is also 
very well matched to K2K physics requirement. 

The tracking efficiency is shown in Fig.|31as a 
function of p and 9 X . The contributions of track 
types I, II and III are normalized individually, so 
that the total efficiency is the sum of the three. 
Notice that the loss in efficiency at low polar an- 
gles for tracks of type I (due to saturation ef- 
fects in NDC1) is largely compensated by tracks 
of type II and III, so that the total tracking effi- 
ciency is rather flat. 

Finally, the yield must be corrected by the pion 
efficiency and purity, which is computed from the 
data themselves, using samples of pions, protons 
electrons and kaons of various energies tagged 
with the help of the beam detectors. The pion 
correction factor defined as the ratio of purity 
over efficiency is shown in Fig. 0] for the three 
track types and for 1.5, 3 and 5 GeV beam par- 
ticles. 

Preliminary studies show an almost diagonal 
migration matrix, where migration is mainly due 
to finite momentum and angular resolutions, but 
not to systematic effects (momentum and angular 
biases, etc). 

Figure shows the pion yields as a function of 
momentum after succesive corrections (e.g., ac- 
ceptance, tracking efficiency and pion correction 
factor). No attempt to deconvolute the data (e.g, 
correct for migrations) have been done yet, and 
no absolute normalization is computed. 
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Figure 2. Acceptance in P for \8 y \ < BOmrad, 
\8 X \ < 200mrad. Acceptance in 9 X for \8 y \ < 
50mrad and P > lGeV 



Figure 3. Tracking efficiency as a function of 
p (left) and 6 X (right). The colors indicate: 
black=total, green=type I, red=type II, blue= 
type III. 



5 




Figure 4. Pion correction factor. The colors 
indicate: red=1.5GeV, green=3GeV, blue=5GeV 



4. Conclusions 

We have presented a preliminary analysis of one 
of HARP data samples, corresponding to a thin 
Al target, at beam energies of 12.9 GeV. This 
analysis is relevant for the K2K experiment and 
a first step towards an upcoming measurement of 
the fully corrected production cross section. 
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Figure 5. Upper; Raw pion yield. Lower; Yield 
corrected by pion correction factor, efficiency and 
acceptance. The red filled histo corresponds to 
the proton+electron background 
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